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1 Introduction and motivation 

The last few years have seen a renewed focus on string dual descriptions of QCD following 
the discovery of the AdS/CFT correspondence [1-3] and subsequent attempts to construct 
holographic string/gravity duals of QCD in both top-down [4, 5] and bottom-up approaches 
[6-10]. These models have been remarkably successful in reproducing some of the low- 
energy features of QCD, but also suffer a number of flaws discussed at length in the 
literature [12-15]. 

For concreteness, let us consider the "hard wall" model of [7]. This model contains 
five-dimensional (5d) fields A'^^ ^, ^ ^^'^ ■^"'^ dual to the operators qi^^f^qL, QrI n't"' Qr 
and q'^q^ respectively. These fields live in a 5d Anti de Sitter [AdS) space with a hard 
cutoff on the radial coordinate at z = zq with I/zq = 346 MeV. The model gives rise 
to a spectrum of 4d fields whose overall scale is determined by the infrared (IR) cutoff 
zq. The lightest modes of definite parity are the vr, p and ai mesons; the model generates 
an excellent fit to the masses, decay constants and coupling constants of these mesons in 
terms of three parameters: zq, the quark mass ruq = 2.3 MeV, and the scale of chiral 
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symmetry breaking a = (308 MeV)'^ (using the Model B values of [7]). Other "bottom-up" 
models [8-10] have a very similar structure. The "top-down" model of [4] is based on a 
D-brane construction in string theory, and whose field content is uniquely as a result. It is 
roughly as successful at predicting the properties of the low-lying mesons as the bottom-up 
approaches. 

Despite its successes, the hard wall model suffers from a variety of issues, some of 
which we address in this paper. The most obvious problem has to do with the structure of 
the excited states in the meson spectrum. The first excited p oi u meson state is predicted 
to lie at m{p') = xq^2/zq = 1910 MeV where j;o,2 is the second zero of the Bessel function 
Jo(x), while in reality the mass of the observed excited state is at roughly 1450 MeV. In 
addition, the mass squared of higher excited states scales with the excitation number as 
v? rather than as n as expected from Regge theory and semi-classical arguments [12, 16]. 

One solution to the problem of the excited meson spectrum, analyzed in [16], is to alter 
the z-dependence of the metric and of the dilaton profile. In [16] a simple modification was 
suggested, which leads to a spectrum of excited vector mesons of the form = S(n + 1) 
where S is a constant related to the QCD string tension. For an appropriate choice of 
S this fits the low-lying spectrum reasonably well and leads to the expected asymptotic 
behavior at large n. 

However, we will see that this cannot be the full story. The hard wall model does 
not even fully describe the light meson spectrum. It gives reasonable predictions for the 
lowest-lying isotriplet vector and axial-vector mesons, the p and ai (and also the isosinglet 
uj and /i if the gauge group is generalized from SU{2)l x SU{2)r to U{2)l x U{2)ji) but 
does not include fields which give rise to the = 1+- hi and hi mesons whose masses, 
at ~ 1200 MeV, lie just below those of the ai and /i. Furthermore, as discussed in [11-13], 
from QCD one actually expects two types of p mesons. The first type couples dominantly 
to the usual vector current qj^f^q, the second to the tensor operator qa^'^t^q. In the 
limit of unbroken chiral symmetry these two operators lie in different representations of 
the chiral symmetry group, leading to two distinct towers of states with distinct (chiral) 
quantum numbers. When chiral symmetry is broken, the operators mix and the physical 
p mesons will be linear combinations of the two p varieties. One cannot expect to get a 
correct spectrum of excited p and to mesons if this structure is ignored. 

A possible solution for including the hi and bi mesons was suggested in [17]; one should 
add 5d fields dual to the interpolating operators for these mesons. These operators are of 
the form qa^'^t"'q; incorporating their dual fields is a very natural generalization of the hard 
wall and other dual models of QCD. Once one includes fields dual to some of the canonical 
dimension three operators in QCD as was done in [7] it is natural to include fields dual to 
all the canonical dimension three operators in QCD. Although the solution is conceptually 
clear, we will need to overcome a number of technical issues in order to implement this 
idea. 

In what follows we will argue that there is a natural extension of the hard wall model 
which includes complex, antisymmetric tensor fields bMN transforming in the bifundamen- 
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tal of the U{Nf)L x U{Nf)R gauge symmetry, that are dual to the operators qa'^'^t"-q^. We 
win show that these fields have a natural first order action, familiar from the description 
of charged tensor fields in AdS^ x [18, 19], and that their expansion in terms of 4d 
fields leads to a tower of mesons with the quantum numbers of the hi and bi as well as an 
additional tower of mesons with the same quantum numbers as the p and uj mesons. We 
will also identify the lowest dimension operator which induces mixing between the two dis- 
tinct towers of p/cj mesons in the presence of chiral symmetry breaking, and which is also 
responsible for decay processes like hi — >• p + ir and bi — )• w + vr. Our treatment follows the 
standard formalism of string/gauge duality. In this we differ from [20] where it is claimed 
that modifications to the formalism involving the asymptotic behavior of fields is required 
in order to incorporate tensor mesons in the presence of chiral symmetry breaking. Our 
description of tensor mesons also seems to be somewhat at odds with the description given 
in [21]. It is not clear to us whether the formalism used in [21] included the possibility of 
fields with a first order Lagrangian. [16] does not address this additional tensor operator, 
and in fact we should note that despite its encouraging fit with the excited p spectrum, 
[16] does not address the hi and hi states. Once we include fields dual to these excitations 
we necessarily generate additional vector meson states, presumably ruining the agreement 
between the excited /o's of [16] and experiment. 

In addition to these specific issues, which are directly addressed in the body of the 
paper, our analysis will also touch on some important problems of principle that arise in 
the construction of QCD duals. These involve the nature the 1/Nc expansion, the validity of 
a local 5d field theory description and the field theory interpretation of the asymptotically 
AdS spacetime. 

Let us consider the last point first. It is often said that the conformal invariance of Anti 
de Sitter space matches the conformal invariance that QCD enjoys in the ultraviolet (UV) 
because of its asymptotic freedom. This is not entirely correct. In the best-understood 
example, the AdS/CFT correspondence, the conformal isometry of AdS space is in ac- 
cordance with the conformal invariance of = 4 SYM at large 't Hooft coupling. It is 
important to note that in this limit A^ = 4 SYM is not a free field theory, but is rather a 
nontrivial CFT with dimensions and correlators which differ from those of free field theory. 
In the analysis of [7] the 5d gauge coupling is determined by matching onto the UV behav- 
ior of the vector current two point function in QCD, which can be computed reliably from 
a one-loop diagram because of asymptotic freedom. This is a very special case, however, 
because the anomalous dimension of the (conserved) vector current vanishes. In contrast, 
the tensor operator qt°'(T^yq is not conserved and has a nonzero anomalous dimension. We 
will see for example that matching its anomalous dimension as derived from the AdS / CFT 
correspondence to the UV behavior of QCD leads to an unphysical value of the mass of 
the lowest bi/hi mesons. This suggests that the correct interpretation of models like [7] is 
that they are modeling QCD as a theory which has a "conformal window" of energy scales 

^We work in the context of the hard wall model for simplicity. A more complete dual description of 
QCD will undoubtedly incorporate the tensor fields discussed here as well as modifications of the metric and 
dilaton profile in order to produce the correct asymptotic behavior of excited states. The dual description 
of the thermodynamic properties of QCD also requires such modifications. See e.g. [22] for a recent review. 
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in which the theory behaves as a nontrivial CFT and then has conformal invariance broken 
at an IR scale l/-2o- This is clearly not the behavior of real world QCD, but apparently for 
some purposes it is a reasonable approximation. With this picture in mind we can obtain 
the appropriate 5d mass of the tensor field (dual to the anomalous dimension of the tensor 
operator) required to fit the observed mass of the bi/hi mesons. 

We turn now to the issue of a local 5d dual description of QCD. In current models one 
studies a 5d action for massless fields such as the gauge fields ^l,_r. In the top-down model 
of [4] these fields arise as the lightest open string modes. One can also include the lightest 
closed string modes: the graviton, dilaton and antisymmetric tensor fields, which provide a 
dual description of the glueball spectrum of QCD. One then either writes down the lowest 
dimension interactions involving these fields, or derives the lowest dimension terms using 
an a' expansion in string theory. In other words, one works in a "supergravity" limit of 
the full string dual. Such an approximation is usually justified only if there is a clear 
separation of scales between the massive string modes and the energy scale being used to 
probe the massless string modes. There is no evidence for such a separation of scales in 
real world QCD. For example, the string scale extracted from the Regge trajectories of 
meson states is a' ~ 0.88 GeV^ which leads us to expect massive string modes at a scale 
of order l/^/a' ~ 1 GeV. Yet in both the top-down and bottom-up models this scale 
is comparable to the scale of excited states of the massless modes. Put another way, in 
the current models one keeps only the lowest dimension operators constructed out of the 
massless modes. If there is no separation of energy scales, there is no reason to expect 
higher dimension operators to be suppressed, or for that matter, no reason to expect that 
we can use a local 5d field theory description at all. 

Nevertheless, in what follows we assume that it makes sense to use a local 5d field 
theory description. One possible theoretical justification for this has been suggested in 
[16, 23]. [23] proposes that there exists a local 5d description in the — )• oo limit for 
small 't Hooft coupling, A = g'^N, even though this is the a' — )■ oo limit rather than the 
a' —7- limit one usually considers to obtain a local 5d action, and that, furthermore, 
the existence of such a local description is linked to the existence of an infinite family of 
conserved tensors of arbitrarily high spin. See [24] for a review. 

In [16] it was argued that one can add 5d fields dual to an infinite family of twist two 
operators in QCD of the form q^(P-^ D^^^ ■ ■ ■ D^'^\ while maintaining general coordinate 
and tensor gauge invariance at quadratic order in the fields. The tensor field bMN that 
we add is also dual to a twist two operator in QCD, and we will show that again one 
can maintain general coordinate and gauge invariance at quadratic order. In fact the full 
leading Regge trajectory of QCD including daughter states seems to be reasonably well- 
described by fields which are dual to the twist two operators of QCD. It is also known from 
semi-classical arguments that operators of fixed twist receive corrections to their anomalous 
dimensions which are smaller than one would naively expect, growing only as the logarithm 
of the spin S [25]. Finally, we remind the reader of one other example where one finds 
remarkably accurate results without the presence of an obvious low-energy limit, namely 
level truncation in open string field theory [26]. 

We will have more to say about higher dimension terms in the effective action when we 



- 4 - 



construct the effective action for tensor fields in section 3. Whether or not these theoretical 
speculations are borne out, our fundamental justification for using a local 5d field theory 
approach is pragmatic. The best way to test any model is to extend it into new regimes, 
perform new calculations, and compare it with data. 

The outline of this paper is as follows. In the second section we review basic facts 
about meson masses and decay constants in QCD as well as the quantum numbers of 
various quark bilinear operators. In the third section we introduce the fields dual to the 
dimension three tensor operators, and construct the free 5d Lagrangian. We then compute 
the spectrum and decay constants in the free theory. In section 4 we discuss interaction 
terms. There are a number of possible Lagrangians that one could consider and the proper 
choice requires a careful discussion of discrete symmetries. We also discuss the role of the 
1/Nc expansion and dimensional analysis. In the final section we conclude and discuss 
some open issues. The appendices contain some useful technical details. In a subsequent 
paper we will present a detailed analysis of chiral symmetry breaking in this model, and 
the effect of the new interaction terms on the spectrum. 

2 The QCD picture 

We begin by describing the features of QCD that are most relevant to our analysis. 
2.1 Operators and two-point functions 

In QCD with two flavors the quark bilinear operators with naive dimension three are 

q{x)eq{x) , 

q{x)f--if,q{x) , 

qix)t''af,^q{x) , (2.1) 

with cr^jy = i{'y^^u — 7i/7/j)/2. The index a = 1,2,3 for SU{2) flavor symmetry or a = 
0,1,2,3 for U{2). In what follows we mostly suppress the flavor indices unless they are 
required for clarity. As the hard wall model of [7] contains fields dual to , O^, , O"^, 
we emphasize here the role of the tensor operator O^. 

The operator is odd under C but contains both parity even and odd parts and 
thus has a non-zero amplitude to create both J^*^ = 1 and 1^ states p^"-\ 6^"^ 

(0|O^"(x)|4")''^(A:)) = ^ft^.^^.e^^'-kf^S'^'^e-^'^- , (2.2) 

(0|Oj^(x)|p(^)''=(A;)) = -!=/J.W(eW^fc, - e^^^k,)6^'^e-^'-^ . (2.3) 

where e^^^s^^^ denote transverse polarizations. In the rest frame of the mesons with 
k = (m, 0), (2.2) implies that the 6i mesons are created by the transverse components Of^ 
while (2.3) implies that the p mesons are created by the longitudinal components O^^. This 



O^'^^ix) = 
O^'^^ix) = 
Ol'^ix) = 
O^'^ix) = 
Ol,\x) = 
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is consistent with the parity assignments of the 61, p. Of course, p states are also created 
by the vector current^ 

(0|O^'«(x)|pW,c(^)) ^ i^in)jV,in)^(n)^^ac^-ik.. (2.4) 

We will see below that there are a number of ways to split up the degrees of freedom in 
, and that different decompositions are practical in different contexts. For instance, we 
can project out the transverse and longitudinal parts of Oj^^ in momentum space using 
the transverse and longitudinal projection operators "P^, (see Appendix A for details) 
- these isolate the hi/hi and a;//3-type states. We define O^^ and O-^"" via ('P"'")/i^O^^ = 

A;20^,^and {vWf.^Olp = k^oll 

On the other hand, when we discuss the transformation properties of under the 
chiral symmetry group U{Nf)L x U{Nf)fi it is convenient to decompose the tensor into 
self-dual and anti-self-dual parts 

Oj:^^ = qcy,.^^q , (2.5) 

which obey 

Oj:i^ = ^^e^^'^Olf . (2.6) 

The operators O-^'^ transform as (1,0) and (0, 1) under Lorentz transformations and obey 
(0'^'+)* = O^-. That is, complex conjugation flips the spacetime chirality, as is familiar 
from the behavior of Weyl spinors. Under the chiral symmetry group 0^'~^ transforms as 
a bifundamental, O"^'^ ~ (N^,Nj) while 0^'~ transforms as the conjugate, (Nj,Nj). 
We can also construct combinations of deflnite parity as 

Oj. = 0++0;, , (2.7) 

which transforms as a tensor under parity and 

0^'" = 0+-0-, (2.8) 

which transforms as a pseudotensor. These are not independent as a consequence of the 
identity <t^i/75 = ie^^Xp'^^'^ which implies that 

= \e,^''Ol . (2.9) 

Our primary source for information about the properties of the dual theory will be the 
two-point function, whose momentum-space poles mark the masses of excitations, while 
the residues at the poles provide the decay constants. The large-Q behavior of holographic 

^The literature contains many different conventions for tlie normalization of decay constants, some of 
which differ from those used here by factors of 2, y/2 and powers of the meson mass. Our normalization for 
the decay constant fp is fairly standard, our normalization for the decay constants /t, and /J agrees with 
that used in [35, 36] once the difference in isospin conventions is taken into account. 
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correlators is often compared directly to QCD. We thus summarize the relevant QCD 
two-point functions here: 



n3?(^) = e*'"(0|O^'^(x)Oj''(0)|0) . (2.10) 

It is conventional to separate out kinematical and group theory factors which are dictated 
by parity, current conservation and Lorentz invariance, so we define 

nS?(^) = (^;.vn^'^"(^^) + vi,,^u^'^Hk')) . (2.11) 

These two point functions receive both perturbative and nonperturbative contributions 
and have been studied in [27-31]. At large —k'^ these can be computed perturbatively in 
QCD and one finds 

n^'^(fc^)^^iog(-A.^), 

n^'^(fe^)^-^m,log(-fc2), 

nT,T||(^2) ^ ^log(-A:2) , U^'^Hk') ^ -^log(-fc') , (2.12) 

where rUg is the quark mass. In the chiral limit, niq — ?■ 0, the tensor and vector currents do 
not mix perturbatively. This is consistent with the observation that 0^,0^ transform in 
different representations of the chiral symmetry group. Of course the perturbative mixing 
is proportional to the explicit chiral symmetry breaking parameter - but one should expect 
nonperturbative mixing, even in the chiral limit, due to the spontaneous chiral symmetry 
breaking of the quark condensate. 

Given (2.2)-(2.4), it is clear that 11^'^ and II^-'^II should feature resonances correspond- 
ing to uj/p exchange, while n^'"^-*- should have resonances corresponding to hi/bi exchange. 
By inserting a complete set of states into the two-point functions one learns 

n^'^(A;)|polcs = -Yl ,2 \ W^2 ' ^^'^ (^)lpoles = Yl ,2 \ (")n2 ' 

n k^ - (ml 'Y 

The decay "constants" here are the same as those appearing in (2.2)-(2.4). In general they 
run with scale. 
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2.2 Low-lying hadron spectrum and decay constants 

As mentioned in the previous section, we can think of these vector and tensor operators as 
generating towers of increasingly massive spin-one resonances from the vacuum. The PDG 
summary table [32] lists three spin-one mesons with the quantum numbers of the isotriplet 
p meson: /9(770), /3(1450), and p{l7m) with masses in MeV of 775.49 ± 0.34, 1465 ± 25, 
and 1720 it 20 respectively; and three spin one mesons with the quantum numbers of the 
isosinglet uj: lo{782), a;(1420), and a;(1650) with masses of 782.65 ±0.12, 1400 - 1450, and 
1670 ± 30. Additional states such as the /9(1900), p(2150) are mentioned in the complete 
particle listings and their mass spacing was quoted in [16] as part of the evidence for 
linearity in n of the mass squared of exited meson states m^, but their existence must be 
regarded as uncertain. In comparison, the experimental evidence for the three lightest p 
states is now quite compelling [33]. 
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92.4 ±0.35 


[32] 
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[32] 


fai 


433 ± 13 


[32] 




1281.8 ±0.6 


[32] 



Table 1. Masses and decay constants of low- lying mesons in MeV. The decay constants fj 
are evaluated at a scale of 2 GeV. Our focus is on 1^^ and 1"' states; the lightest axial- vectors 
and pseudoscalars are included for completeness. 

We will denote the lowest three mass eigenstates by p, p' and p" . The corresponding 
vector and tensor decay constants, as defined by (2.2)-(2.4), will be denoted fp, fp', fp" and 
fj , fj,, fj/,- In table 1 we summarize experimental results for masses and lattice results for 
decay constants of the tensor and vector mesons. 
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3 Holographic dual at the free level 



Having reviewed the experimental and lattice data on the QCD side, we now extract this 
data from a holographic model. We first identify the field dual to the tensor operator 
and its Lagrangian on the hard wall background. We then compute the masses and decay 
constants from the free equations of motion. 

3.1 Field-operator correspondence and quadratic action 

The gravitational background of the hard wall model consists of AdS^ truncated at a finite 
radius [7]. Working with coordinates x^'^ = {x'^,z) we use the metric 

f 

ds'^ = ^(r].,ydx>'dx'' - dz"^) , e < z < zq , (3.1) 

where i is the AdS radius and zq the infrared (IR) cutoff. (The AdS radius is often set 
to one, but we find it both useful and less confusing to keep it explicit). We will perform 
calculations at finite e, taking the limit e — )■ at the end. 

In what follows, q^^^pQ^ will denote the 5d Levi-Civita tensor with normalization 
gOi232 _ i/y/^. "VVe will often rewrite parts of the 5d action in terms of 4d fields; when we 
do so Greek indices /x, i/, • • • will be raised and lowered with the flat Minkowski metric rj^y. 
We also use e^^'^^P for the 4d Levi-Civita symbol with e0i23 _ 

The hard wall model of [7] introduces the fields dual to the dimension three operators 
O^'^ and to O^'^ in the above background. These are, respectively, a tachyon field X 
in the bifundamental of the U{Nf)L x U{Nf)fi gauge group, and the (axial)-vector gauge 
fields A, V. The fields live in a background geometry (3.1). 

We now extend the hard wall model to include a field dual to the tensor operator . 
The obvious choice is a two- form potential bMN- To match the global flavor transformation 
properties of O"^, b^jN must transform as a bifundamental under the U{Nf)L x U{Nf)ji 
gauge group, and must therefore be a complex field. 

In order to match the six physical degrees of freedom generated by (3 for each of 
the two massive vector-like states), bMN should also carry a total of six degrees of freedom 
on shell. The complex field bMN has twenty independent real components, but as will be 
seen explicitly later, the equations of motion imply a Proca-like condition which removes 
eight of these, leaving twelve independent on-shell degrees of freedom - twice as many as we 
want. In order to eliminate these additional states, one might consider imposing a tensor 
gauge invariance, but this turns out to be impossible. Since b transforms nontrivially under 
the gauge group, the usual three-form field strength H = db is not gauge covariant. On 
the other hand, the gauge covariant antisymmetrized derivative 

'HpMN = '^D[pbMN] (3-2) 

is not invariant under tensor gauge transformations 6b mn = "^^IM^N] and as a result one 
cannot use tensor gauge invariance to remove the unphysical degrees of freedom. 

This problem is familiar from the AdS/CFT literature where one also finds two-form 
tensor fields transforming under the 50(6) gauge group in the dimensional reduction of 
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IIB string theory on [37, 38]. These fields are described by a first-order Lagrangian 
discussed in [39] and have been analyzed in the context of the AdS/CFT correspondence 
in [18, 19]'^. Though there is no tensor gauge invariance, none is needed since a first-order 
action has fewer independent degrees of freedom. In a first-order formalism half of these 
are momenta and half are coordinates, so one is left with six independent real degrees of 
freedom. This correctly matches the six physical degrees of freedom (two massive spin one 
vectors) created by the dual operator . 

According to the AdS/CFT map, a (i/2-form field in AdSd+i obeying a first-order 
equation with mass /i (in units of the AdS radius) is dual to an operator with dimension 

A = i(d + 2|/i|) . (3.3) 

In contrast to the scaling dimension of the conserved current Jy, we should expect correc- 
tions to the naive value of A = 3 for . We thus leave fi arbitrary. 

We now have the following extended version of the hard wall action: 

S = 5hw + 5*08 + Ssd + Sint , (3.4) 

where 

5hw = I d'x^ tr + - ^(F| + f2)| (3.5) 

is the action used in [7], and 

5cs = ^ ^ (u;,{Al) - u.iAnij , (3.6) 

with trF^ = duj5, is the Chern-Simons term needed to match the flavor anomalies of QCD 
(see e.g. [41]), and 

is the first order action for the antisymmetric tensor field written in terms of differential 
forms^. We use a bar to denote Hermitian conjugation on group indices-e.g. X = . Since 
we work in a basis of Hermitian generators, this amounts to complex conjugation of the 
individual flavor components. The interaction term Sjnt will be discussed later. Note that 
we have explicitly inserted factors of the AdS^ radius, £, in such a way that the couplings 55 
and gb are dimensionless. The mass dimensions of the flelds are [X] = [b] = 3/2, [^l,_r] = 1. 

First order actions on manifolds with boundary require the addition of boundary terms 
as dictated by the requirement of a consistent variational principle [42] or from consistency 

^At the level of free field theory one can establish an equivalence between the first order formalism and 
a second order formalism with a conventional kinetic energy term and a Chern-Simons mass term as in [40] 
but is is not clear that this equivalence can be established when interactions are included. 

■'We have b = ^buNdx^' dx^ and Db = ^Dypbmi^^dx^ dx'^^ dx^ , and the Hodge dual is {*b)MNP ~ 
l^MNp'^^^QR- Since 6 is a bifundamental, like X, the gauge covariant derivative acts as DpbmN ~ 
dphniN — iAL,pbMN + ibMN Ar^p. 
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when passing between the Lagrangian and Hamiltonian formulations of the theory [43]. 
These terms play a crucial role in giving a precise definition to the AdS/CFT correspon- 
dence since the bulk action vanishes on the equations of motion, leaving the boundary term 
to generate correlation functions when we vary with respect to the sources. 

The boundary term in (3.7) was first found by [18], but we can also obtain it via a 
variational argument as in [42]. Since this term will play an important role, we review 
the argument here. The boundary term is required for the variational problem to be well- 
defined: without the boundary term the variation of the action evaluated on shell will 
not vanish, implying that solutions to the equations of motion do not represent stationary 
points of the action, and thus invalidating the stationary phase approximation to the 
partition function. 

We begin by analyzing the equations of motion to determine the appropriate degrees 
of freedom. Taking the variation of Sgd with respect to b, the (free) equations of motion 
are 

(d-i^*)6 = 0. (3.8) 

We review the solution to these equations in some detail in Appendix B. Here we note 
simply that b^z is determined in terms of 6^,^, and one may derive separate equations for 
the self-dual and anti-self-dual parts of 6^,^. Writing 

bf,u{x, z) = 6+^(x, z) + b~^{x, z) with b^^ = ±^e^yp„b^P'' , (3.9) 
the (anti) self-dual pieces obey the equations 

6±(fe,z) = 0, (3.10) 

where b{k,z) = f d!^xe'^^'^b{x, z), as usual. Near the UV boundary z ~ e — )• 0, solutions to 
(3.10) behave as 

b^u ~ 'S^^.e-^' - , (3.11) 

b~^ ~ A^,^e^~^' - ~a^,^e^ , (3.12) 

where S", s and A, a are self-dual and anti-self-dual polarizations, respectively. Since the 
equation of motion (3.8) is first order, however, these coefficients are not independent. One 
may derive the relation 

(5,s)„^. (3.13) 

(i,a)a/3, (3.14) 
which is consistent with (3.13) since 

(p-L _ y>\\^2 ^ ;^4j^^ Observe that when /i > 0, 5* 
encodes the leading behavior of b^^ near the UV boundary, while when < 0, a encodes 
the leading behavior. It will be important to distinguish these two cases below. 



1 



M/^±2) 



Equivalently, we have 
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Now consider the variation of the bulk part of Sgd, focusing on terms with z-derivatives: 



i 



'ski I {Wu^bpa + bjTuSb;, - b^Jhp^ - b-j'^y^ . (3.15) 

Evaluating the variation on shell, the bulk terms vanish by the equations of motion. The 
relation (3.13) implies that we cannot simultaneously fix 6+^ and at either the UV 
boundary or the IR boundary-to do so would overconstrain the system. At each boundary 
we may fix only one or the other. The natural choice at the UV boundary, from the 
AdS/CFT point of view, is to fix 6+ (fe~) for > (/U < 0). (In the language of [18], for 
fi > 0, b^ plays the role of coordinate and b~ that of canonical momentum, while for /x < 
their roles are reversed). Either way, half of the boundary terms in the last line of (3.15) 
will remain. If we constrain 6"*", say, then the 6b~ terms will be nonzero. 

There is a way out of this quandary. First let us suppose /x > 0. Notice that the 
on-shell value of (3.15) may be written as 

+ / d'xe>^'"'''^T - b^uSb;^y° , (3.16) 

or, moving the first term to the left and using the (anti) self-duality of b^, 

^ (^''^"'' 'wjau'" ^'^'''^'^'^0 {^'^'^'^^ - ^-^^p-ir • ^^-^^^ 

If /i < we should instead isolate 5b~ . We find that the on-shell value of (3.15) may also 
be rewritten as 

^ (^'d'" ^wJ'dM^' 1^'^'^^'^}) = / d'xe^"'"^^ {Ku^'^^ - Wu5h~^Y'^ ■ (3.18) 

Thus, by adding the boundary term, as has been done in (3.7), we allow for a consistent 
variational principle. The on-shell variation of the full (bulk + boundary) action will vanish 
if the boundary value of 5"'" {b~ ) is held fixed in the case /i > (// < 0) . 

We in fact have two choices for the boundary condition. Consider the /x > case for 
concreteness; analogous comments apply in the /z < case. When > 0, the on-shell 
variation of (3.7) will vanish if we either hold 6"*" fixed or set 6~ = on the boundary. 
We will refer to these as the Dirichlet-type and Neumann-type conditions, respectively, 
drawing on the coordinate and momentum characterization of b^ given in [18]. On the UV 
boundary AdS/CFT dictates that we take the Dirichlet-type condition, holding b'^ fixed, 
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since 6+ represents the source for the dual operator. On the IR boundary, the correct 
choice is not immediately obvious. 

The question is compounded by the potential presence of higher-order terms localized 
on the IR boundary and not included in (3.4). In general we should expect such terms to 
appear as a result of "integrating out" nontrivial IR dynamics up to the scale of the IR 
cutoff 1/^0 . These terms could potentially modify the Neumann- type condition to a more 
general mixed condition, which may even be nonlinear, with parameters depending on the 
details of the localized terms^. 

These issues were dealt with in two different ways in the original hard wall model [7], 
depending on the dual field in question. In the case of the scalar field, X, all of these 
unknowns were packaged into a single parameter, the quark condensate, which was taken 
as an input that could be fit to data. In the case of the vector gauge field, V, the simple 
Neumann condition, ^^^(zo) ~ ^^V^l^o = 0, was chosen, which can be motivated as follows. 
First, one can definitively choose in favor of the Neumann condition over the Dirichlet one, 
^V\zo = 0, by requiring that the boundary condition be gauge invariant. This does not rule 
out the possibility of higher order terms in the field strength localized at the IR boundary, 
but the leading term of this sort was considered, and its effect was found to be small in 
practice. In general, one expects terms involving higher powers of the gauge field-on the 
boundary or in the bulk-to be suppressed by powers of 1/Nc, while for X this is not true. 
We will review these arguments in section 4. 

Returning to the case at hand, one may also use Nc counting to argue that higher order 
terms in b are suppressed. We do not have an analog of the gauge principle, however, to 
help us decide between the Neumann or Dirichlet-like condition. In the next subsection we 
use the holographic model (3.4) to compute the tensor-tensor two-point function, (O^O^), 
considering both IR boundary conditions. We will show that there is a strong physical 
argument for choosing the Neumann-like condition over the Dirichlet one: the Dirichlet- 
like condition leads to a massless divergence in the tensor-tensor two-point function, and 
there are no massless particles with the quantum numbers of the /ii/6i-mesons in QCD! 
The Neumann-like condition, on the other hand, gives a sensible finite result for {0^0^) 
in the A: — )■ limit. 

3.2 On-shell action and tensor-tensor two-point function 

Let us now use the quadratic action (3.4), and in particular (3.7), to study the correlation 
functions and the normalizable spectrum of the complex two-form b^jN- The two-form bMN 
corresponds to a total of six real degrees of freedom. We can package these in a variety of 
ways: in terms of real and imaginary, longitudinal and transverse, or self-dual and anti- 
dual parts of b^i,^. It will therefore be convenient to work as much as possible in terms 

^We take the point of view that effects of higher order boundary terms modify the boundary conditions, 
such that a well defined variational principle is maintained. An alternative approach is to hold the simple 
Neumann-type boundary condition fixed, so the presence of higher-order terms leads to a contribution to 
the on-shell value of the action from the IR boundary. On a practical level one may take either point of 
view; these are just two different ways of encoding the effects of such higher-order terms. 

®The components can be eliminated through their equations of motion. 
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of projections onto these parts, and T^H'"*", defined in Appendix A. The intermediate 
steps of the analysis differ between the cases fi > and // < 0, though the final result for 
the tensor-tensor two-point function can be given in a simple form, valid for both cases. 
In order to streamline the discussion in this subsection, we will present the /i > case 
only, and at the end quote the result for general /x. Details of the analysis are relegated to 
appendices B and C. 

We first determine the bulk-to-boundary propagator of b, which we need in order to 
compute two-point functions of the dual operators. 

For the free two-form in AdS^ it is most convenient to work in terms of the self-dual 
piece which satisfies 

b^u = y^u'^b^^ ■ (3.19) 

As discussed by [18], it is the self-dual piece which sources the dual operator: for the 
value /i = 1, which should give the scaling dimension in a conformal theory (i.e. on a 
background that is simply AdS^), it indeed has the appropriate near-boundary scaling 
behavior to correspond to a tensor operator. In the confining theory, of course, the value 
of n should receive quantum corrections. 

The general solution to the free equations of motion for b is reviewed in Appendix B. 
In momentum space one has b^u = b^^^^ + with 

b'l;^{k,z) = S^u{k)zJ-f,-i{kz) + s^y{k)zJ^+i{kz) , 

b~^{k,z) = Af,,,{k)zJ^^+i{kz) + d^u{k)zJ^^i{kz) , (3.20) 

where the anti-self-dual polarizations are related to the self-dual ones by 

■^pu — ~ ^ )fj.u^a![^ ~ ^t^f ~ ykpk^ S pi, — kjjk^ S pp^j , (3.21) 

and similarly for d in terms of s. Meanwhile, 6^^, which plays the role of Lagrange multi- 
plier, is given by bpz = —j^^^fi'^^^k^bpa- The solution (3.20) is appropriate in the generic 
case of non-integer fi; if fi is an integer then the Bessel functions J~p^i should be replaced 
by Yp±i. 

We are interested specifically in the bulk-to-boundary propagator: that is, the solution 
to the equations of motion with the boundary condition that the self-dual field approaches 
a self-dual source on the UV boundary. 

We first impose an IR boundary condition aX z = zq to fix s in terms of S. As we 
discussed above, there are two possibilities for a consistent variational principle: we may 
choose the Dirichlet-like condition where we hold 5b'^ fixed, or the Neumann-like condition 
where we set b~j^{zQ) = 0. Holding 6+ fixed in this context means setting 6^jy(2o) = 0, since 
there is no natural constant antisymmetric two-tensor living on the IR boundary. After 
imposing one of these, we find that the solution takes the form 

b^^{k,z) = Spy{k) [zJ^p_i{kz) - {k, ZQ)zJi^,+i{kz)\ = S^y{k)B^{k, z) , 

b'l.uik^z) = Apy{k) [zJ-i,+i{kz) - c'^{k,zo)zJf,-i{kz)] = A^^{k)B~{k,z) , (3.22) 
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where 



(3.23) 

■/-M+i(fczo) A-^zn^ - n 



(The ">" labels are a reminder that these expressions are appropriate for the case ^ > 0). 

The UV boundary condition fixes S^y in terms of the source for the dual operator. As 
we approach the UV boundary, 2 = e — )■ 0, the leading scaling behavior of 6 is ft"*" oc e~^. 
We define the self-dual source, S^y{k), through 



m-l/2 

h+,{k,e) = ^^S,,{k) . (3.24) 



The factors of the AdS radius i have been inserted on dimensional grounds. 6 is a dimension 
3/2 field, while S should have naive dimension 1 since it sources a naive dimension 3 
operator. Using this boundary condition to eliminate 5 in favor of the source S, we arrive 
at our final expression for the bulk to boundary propagator: 



B>{k,z) 
eA'S+(/c,e) 



(3.25) 



Pulling out the explicit factor of as we have done in (3.24), ensures that our final 
expressions for two-point functions will be e-independent and corresponds to working with 
"rescaled" operators. Another procedure commonly appearing in the literature is to simply 
set i^/%+{k,e) = S{k), in which case our result for the two-point function would include 
an overall factor of e^'^, indicating its scaling behavior as one approaches the UV boundary. 
In the e — 7- limit one should trade this factor for a renormalization scale Mr- In (3.24) 
we have chosen to identify the renormalization scale with the AdS radius, Mr ~ This 
is natural since the dimensionful couplings i^^'^g5, i^^'^gt appearing in (3.4) were already 
(implicitly) defined at this scale. 

To summarize, three different length scales appear in the model: the locations of the 
UV and IR boundaries at e and zq, respectively, and i the AdS radius. In terms of the 
dual field theory, the mass scale is the UV cutoff we would include in the computation 
of bare n-point functions. Adding the appropriate counterterms to the action is essentially 
equivalent to replacing with £~^, to give a finite (renormalized) result. is the 
renormalization scale, while Zq^ is Aqcd- 

These statements may seem trivial when we are working with a simplified model con- 
sisting of AdS^ with cutoff. While this model captures basic features such as confinement, 
couplings and dimensions of operators do not run with scale. In a more realistic model this 
running would be encoded in a nontrivial z-dependent geometry. Physically, the picture 
is that the asymptotic part of our AdS slice is not describing asymptotically free QCD, 
but rather is providing an approximate description in a window of scales where the QCD 
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coupling is finite, but running slowly with scale. We have already stressed that the tensor 
operator is not a conserved current, and at strong coupling one should expect 0(1) 
corrections to its charge and dimension - or to /i in the dual language. The renormal- 
ization scale represents a typical scale in this window. The ratio of this scale to the 
QCD scale, parameterized by the ratio is a dimensionless parameter in our model, 

and we will see that it naturally appears in physical quantities like decay constants. 

To compute correlators according to the usual AdS/CFT description, we evaluate the 
action on the solution (3.25), functionally differentiate with respect to the source, and then 
take the e — )• limit. S^^ couples to the self-dual operator, 0^1 in (2.5). It is related to 
the source, T^u, for the tensor operator via 

V = V + ^V'"'V ^ %u{k) = ]^{S^,{k) + S^,{-k)) . (3.26) 

The sole contribution to the on-shell action comes from the UV boundary term, which we 
find to be 

This action leads to the following matrix elements. We have that n'^'-^H = —11-^''^^, as 
is evident from (3.27), and 



n^'^^(A:) 



- 22.-2g2'^(^) Cfc>(fc, zo){k^f^'-^ , /X > non-integer 



96 i 

'-^^ [nc>{k, zo) - log [kH^)] {Uf^-^ , ^ > integer , 



(3.28) 



where c^{k, zq) is given by (3.23) in the non-integer /i case and by (3.23) with J-^^i — )• y^±i 
in the integer /i case. Some of the intermediate steps involved in obtaining (3.27), (3.28), 
are presented in Appendix C. 

Let us consider the extreme IR limit of this result, A; — t- 0. The tensor-tensor correlator 
should be finite in this limit, since there would be no particle interpretation in QCD for 
such a massless divergence. The leading k"^ behavior of c^(/c,zo) is (for both integer and 
non-integer ^) 

hmc,(^,.o) = |^(^_2,+2)^ 6-(.o) = 0. ('-'^^ 

We conclude that the Dirichlet-like boundary condition, 6^(zo) = 0, is physically unaccept- 
able since it leads to a divergence in 11^'"^ at large distances. The Neumann-like condition, 
6~(zo) = 0, on the other hand gives a physically reasonable result for the /c — t- limit of 
n^''^. Per our discussion in the previous subsection, we take the Neumann-like boundary 
condition to be our IR condition for the 6-sector. 
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The final result for the tensor-tensor two-point function, in the case > 0, computed 
using the free dual (3.7), is 



^ > non-integer 
, i-i > integer . 



(3.30) 



In appendices B and C we also analyze the /i < case. The final result is identical to 
(3.30) with fi — )• up to a sign whose origin may be traced to the sign on the boundary 
term in (3.7). Thus the general result, valid for any 7^ 0, is 



n^.5"^(A;) = { 



-22|M|-2g2r(|^|)l«'«J J|„|_l(fc20) 



02 1 .1 2T^m!i I .M (fc^)''^'"' vr ^'-i-\^'^°i - log (kH^) 

^, 22|M|-2g2|^||(|^|„]^)| V ; 1^ J|^|_i{A:2:o) & { 

Next we turn to a discussion of the physics contained in (3.31). 
3.3 Interpretation of results 

The original hard- wall model depends on the parameters zq, niq, and a, as well as 55 
which is expressed in units of the AdS radius i. Our posited extension of the hard-wall 
model contains new parameters fi, as well as explicit dependence on £. All physical 
quantities we compute will be given in terms of these parameters, so in order to make 
additional predictions, we need to fix them using real QCD results. Of course it is best, 
when possible, to fix these parameters exactly using explicit computations in QCD. For 
example, quantities which are not renormalized, such as the anomalous dimension of a 
conserved current, give the same result at weak or strong coupling. The authors of [7] 
successfully employed this strategy to find (75 in (3.4), by comparing the gravity dual and 
perturbative QCD results for the large Jy-Jv correlator. Naively, one might want to fix 
fi and gh by comparing the UV (large momentum) limit of the tensor two-point functions 
to large-momentum QCD. While this method is justified for finding g^ using the two-point 
function of the conserved current Jy, the operator is not conserved, so we have no 
reason to expect that fi and gi, will not receive significant corrections at strong coupling. 
(Incidentally, the issue applies to the scaling dimension of qq, dual to the 5d mass of the 
tachyon field X which was held at its naive value in [7]). 

The two-point function (3.31) encapsulates much of the new information we obtain 
from adding the non-interacting two-form field: the locations of its fc-space poles mark the 
masses of mesonic resonances generated by while the residues at the poles give the 
corresponding decay constants. We can expand the two-point function around its poles, 
and by comparison to (2.13), read off these quantities''. 



non-integer 
, /X integer . 



(3.31) 



^This is equivalent to identifying the masses as the eigenvalues of normalizable modes and the decay 
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By inspection of (3.31) we see that the two-point function (for both integer and non- 
integer /i) has an infinite set of simple poles for values oik = m such that J\^\-i{'mzQ) = 0: 



these define the masses of the states created by O to be m„, = x 



At|-l,n^O 



where x 



\^l\-l,n 



denotes the n-th zero of J\^\-i- Note that each pole corresponds to two degenerate states: 
a 1"' (/ii/6i-like) state and a 1 (w/p-like) state. 
Taylor-expanding near the first pole, we find 



n 



sgn(/^) 



220 

220 



-l,l'^{\^^\)J'\n.\-li^\^l,\--^,l) k'^-mf 



+ 



fi non-integer, 



271 Yl 



|^i|!(|At|-l)!x|,j|_i_i J|'^l_j(x|^l_i_i) fc2-m2 



+ 



integer, 



(3.32) 



where J' is the derivative of the Bessel function with respect to its argument. On comparing 
this with (2.13), we learn that /i should be taken negative, and that the decay constant for 
the first resonances has the form 



gb \za 



9b V^O 



H-1 



H-1 



2x1 



ir(-|/i|)Ji_| |(a: 



1-ImI 



z, 



-1 







2xi, 



-l,l^^Y\^L\-l{^\^i\-l,l) 



m|-i,i) 

/J, non- integer, 
1 



(3.33) 







/Li integer. 



We emphasize that fji^\k) depends on the momentum, and that the prediction made here 
is for the on-shell value of the first resonance. Note also that the sign of /i turns out to be 
significant: we must choose one sign over the other to have real decay constants^. 

Based on the non-interacting model we have described so far, we were able to make 
concrete predictions for the masses and decay constants of the /p' and the hi/bi meson in 
terms of the parameters fi, g^, and zq/I. Note that if we set |//| = 1, its value in perturbative 
QCD, we would find that both the uj' / p' and the hi/bi states are degenerate with the u/p 
mesons generated by the vector current, whose experimental mass is approximately half of 
the Lo' / p' mass! 

Conversely, we can use the measured value of the hi/bi mass of m = 1230 MeV to 
estimate that \p\ runs to ~ 1.82 at this scale, with 

0.82 



1 



9b \zo 



2092 MeV . 



(3.34) 



constants as derivatives of the eigenfunctions. Especially when one studies interactions it is often useful to 
work with the normalizable eigenfunctions, and to write down an effective 4d action for the corresponding 
resonances by integrating out the holographic z direction. We consider this point of view in Appendix D, 
but here it is equally convenient to just work with the two-point function. 

*The fact that fi turns out to be negative - which implies that we are fixing the value of an anti-self- 
dual source on the UV boundary - has no physical meaning, however: if we had instead chosen b as our 
fundamental field and performed the entire analysis above for b instead of b, we would have instead found 
that we had to fix a self-dual source on the UV boundary. 
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While we could use this data to now fix the value of another parameter, we see that 
our description of the tensor operator is still incomplete. So far we have introduced 
with a non-interacting Lagrangian, reading off the masses and decay constants of the 
resonances it creates; these particles are totally ignorant of chiral symmetry breaking, and 
furthermore have no decay modes. In order to make better contact with QCD, lifting the 
degeneracy between the 1^ and 1 states and introducing interactions that allow the 
vector- and tensor-generated u/ p mesons to mix, we must add 5d interaction terms to the 
bulk Lagrangian. We will now discuss how to identify these terms. 

4 Interactions: discrete symmetries and Nc counting 

In order to classify the modes of hMN corresponding to various meson states, we must 
carefully understand the holographic equivalents of discrete symmetries in QCD. Invariance 
of the gravity dual action under these symmetries will also constrain the interactions we 
can legally include in S'int- 

Consider the pure gauge action, ignoring the Chern-Simons terms: 



W = j d'x^til^Fl + F|) . (4.1) 



This sector of the theory has five Z2 symmetries. One of these is T, time reversal invariance, 
which will not be needed in what follows. There is also five-dimensional parity, which acts 
to reverse the sign of an odd number of the spatial coordinates; we take its action to be 

P5 : A^\x^, X, z) ^ (-l)^(^'^) A^'^(x°, -X, z) , (4.2) 

where I{M) is if M = 0,z and 1 if il/ = 1, 2, 3. In addition there are three Z2 symmetries 
which do not act on the coordinates: 

Cl-.Al^ -Al , 

Cr-.Ar^ -A*r . (4.3) 
Regarding the last two, we write the Lie algebra of the gauge group as 

[r",T^] = i/«^^r% (4.4) 

with the T"" Hermitian, normalized to tr(T"r^) = S"'^/2 and the structure constants /"'"^ 
real. It is then easy to check that — (T*^)* obey the same Lie algebra as the so the 
transformation T'^ — >• — (T")* is an automorphism of the Lie algebra. For SU{N) this 
transformation takes the generators in the N to those in the N. For U{1) it just takes 
A^ —7- —Afj^ which is the usual action of C on the vector potential, so it makes sense 
to call such a transformation charge conjugation. The action S'gaugo has separate charge 
conjugation symmetries for U{Nf)i and U{Nf)R. 
We now consider adding the Chern-Simons term 

5cs = ^ / {^,{Al) - o^5(^r)1 . (4.5) 



M 
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This term changes sign under P5 which changes the orientation of M and it also clearly 
changes sign under P. Thus this term is only invariant under the combination P = P5P. 
This is analogous to the discussion in [44] of Z2 symmetries in the pion low-energy effective 
action where one finds an extra Z2 symmetry of the action which is not a symmetry of 
QCD, and which is only broken by the addition of the Wess-Zumino-Witten term. 

Since trF| = du!5{AL) and under Cl we have trF| — )• —iiFj^F'£F'£ = — trF| we see 
that uj^{Ai), uj^{Api) are odd under C^, Cr. Thus the Chern-Simons term is invariant 
under the combinations ClCrP, or under ClCrP. We will pick the first one and call it C, 
C = ClCrP. Thus the pure gauge action, including the Chern-Simons term, is invariant 
under a x Z2 symmetry generated by C and P. 

We now extend these symmetries to the bifundamental fields X, h. The covariant 
derivative term is invariant provided that C : X — )• ibX^,P : X — )• itX^ = ±X, 

but since changing the sign of X can be accomplished by a gauge transformation we can 
choose the action to be 

P : X ^X . (4.6) 

where the action of P on the argument of X is the same as P^ in (4.2) and will be suppressed 
from now on. Note that if we follow the treatment in [7] and expand X around its vacuum 
expectation value as X = Xq{z) exp(2i7r"r") then the above action of C, P agrees with the 
canonical assignments of C,P to the pion fields^. The terms in the action involving b are 
invariant under C, P provided that 

C : buN i^MAT ) 

P : bMN ^ ±(-1)^(*'^'^)6m7V . (4.7) 

where I{M,N) is the number of indices in (MN). We can fix the sign in the action 
of C by comparison to the standard assignment of C = — 1 to the hi meson. The action 
of P on the real and imaginary parts of bMN differs by a sign, so the choice of sign in the 
action of P just changes which fields we associate to the real and imaginary parts of bj^N- 
In what follows we choose the plus sign so that 

C : bMN —bJiN ) 

P : bMN ^ (-l)'(^'^) W . (4.8) 

Any bottom-up 5d dual contains, in principle, an infinite number of possible interaction 
terms. Such terms are not only relevant to higher point functions between QCD operators: 
when qq acquires a vev, breaking chiral symmetry, these interactions also contribute to 
two-point functions, modifying the location of the poles and values of the decay constants. 
In the original hard wall model [7], this mechanism breaks the degeneracy between vector 
and axial-vector modes. Once we include bMN, such terms will induce the mixing between 
vector meson states generated by , and those generated by . 

^Later on we will find it more convenient to make a field redefinition which embeds the 4d pion field 
entirely in the axial- vector gauge field. 



-20- 



It is not trivial to decide which of this infinite set of terms to include in our analysis. In 
the language of effective field theory we often group interactions by dimension and assume 
that higher dimension operators are suppressed by some scale. From the action in (3.4) we 
find that the mass dimensions of the fields are 

[X]=3/2, [6] =3/2, [F]=2. (4.9) 

so in principle we should only consider interaction terms, such as those roughly of the form 
bFX, which have dimension 5. No real separation of scales exists in the hard wall model, 
however, so it is not clear that this classification of dominant interaction terms is valid. 

Another possible organizing principle is to consider instead the large- A'c scaling of 
interaction terms, and to include only those leading in 1/Nc- We first recall the standard 
large Nc counting rules for QCD as reviewed in [45, 46]. Consider quark bilinear operators. 
Hi, which can include any number of gauge fields as well. The large Nc rules for QCD 
diagrams give 

{Hi---Hr)^N'c^^ . (4.10) 

The operators Hi are normalized so that ^/NcHi has an amplitude of order one to create 
a one meson state. Therefore matrix elements of operators that create mesons with unit 
amplitude behave as 

(v^^i • • • ^/WcHr) ~ iV,^-"/2 . (4.11) 

The two-point function is 0(1) and the three-point function is 0{1/ ^/Nc). This shows that 
the coupling is 0{\/y/Wc), and a two-body decay amplitude for a meson is 0{l/^/N'c) so 
mesons become stable in the large Nc limit. 

In the hard wall model we have been using a Lagrangian schematically of the form 



S ^ }^ f F^ + ^ I I {\DX\^ + |X|2) + . . . , (4.12) 
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where • • • represents additional interaction terms. If we define rescaled gauge fields A' via 
A = g^A' and a rescaled tensor field via h = g^h' then we have 



S ~ j F + J b" + J {\DX\' + \X\') + ■■■ . (4.13) 

By computing two- and three-point functions and matching to the Nc counting of the 
field theory, one finds g^ ~ \/\/Wc which is consistent with the result found in [7] from 
matching the UV behavior of the vector current two-point function. The field b only appears 
quadratically in what we have done so far, so the Nc dependence of g^ is not determined. 
If we use the above Lagrangian to compute two-point functions by the usual prescription 
in AdS/CFT they will be order A*"^ in agreement with (4.11). Three-point functions, which 
involve couplings between three gauge fields, or a gauge field and two pion fields and so 
on, will involve a factor of 55 from the rescaling, and so will indeed be of order 1 /^/Nc- 

Terms involving the tachyon field X are more subtle. Consider for example a coupling 
of the form bXF which can contribute to two-point functions when X is set equal to 
its vev, and apparently also to three-point couplings when we include fluctuations of X 
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(pion modes). If we only focus on pion couplings then we can write X = Xge^*'^ which 
we leave invariant while redefining the gauge fields and h field by what would be a gauge 
transformation if we were transforming X as well and with Vl = = e*'^- This removes 
the pion field from X and puts it entirely in the axial gauge field. After this redefinition, 
we can replace any X appearing in the action with Xq. Now the bXF term contributes 
only to the two-point function, and one might conclude that its coupling should be of 
order N^. However in general one could also look at fluctuations in the magnitude of X 
which are dual to the broad a resonance of QCD. This term then contributes to three-point 
couplings involving cr, a vector meson and a /11/61 meson and one concludes that in fact 
its coupling should be of order 1/ ^/Wl. This term thus gives a subleading contribution to 
the two-point function in the 1/Nc expansion. Terms containing even more powers of X 
will be suppressed by additional powers of l/y/Wl and should be small at large Nc- 

We now consider two interaction terms that are consistent with all global and local 
symmetries and are the leading terms that lead to tensor- vector mixing in the uj/ p sector, 
and break the degeneracy between the and tensor-w/p spectrum. While conceptually 

straightforward, a thorough analysis of the interacting model turns out to be technically 
complex, and is postponed to a forthcoming paper [47]. 

The first term we consider is the unique cubic, dimension five operator which is P, C 
and U{Nf)L x U{Nf)R invariant: 

Sg, =gij d^'x^tv {bMNF^'^'X + XFl'^^hMN + XF^'^Imn + HinFI'^'x} . (4.14) 

— 1/2 

By the earlier argument after rescaling fields we find gighdb ~ 0{Nc ) which implies 
9i9b ~ 0{N^). Evaluating X on its chiral symmetry breaking vev, {X) = {X) oc 1^(2;), 
one finds a quadratic order term that mixes the vector gauge field V with the two-form b. 
The mixing modifies the mass spectrum of both vector and tensor-oj//) states, and breaks 
the degeneracy between the spectrum of hi/bi and tensor-w/p states. 
The second interaction term we consider is the dimension six term 

Sg, = g2i j d^x^ti {bMNXb'''''X + bMNXb'''''X] , (4.15) 

— 1/2 

and we have 5256 ~ 0{Nc )■ Although this term does not contribute directly to tensor- 
vector mixing in the oj/ p sector, it does contain quadratic terms, when X is evaluated on 
its vev, that break the degeneracy between the hi/bi and tensor-w//? sectors. 

In a follow-up paper we will analyze the interacting model (3.4), with 5int = Sg^ + Sg^. 

5 Conclusions 

We have presented a consistent formalism for including fields dual to the antisymmetric 
tensor quark bilinear in a five-dimensional gravity dual description of QCD. This op- 
erator creates J^^ = 1+- mesons such as the hi and bi mesons as well as J^^ = 1 
a;/p-like mesons. 

While the principle of including this new field in the holographic framework is rela- 
tively straightforward, the implementation required several novel strategies. One was the 
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use of a first order action, required for gauge invariance under U{Nf)L x U{Nf)^ flavor 
transformations while demanding that the fields give rise to the correct number of degrees 
of freedom in four dimensions. Another key technical point was to determine a unique IR 
boundary condition for the new field, which was based on the physical requirement that 

should not generate a zero-momentum mode. We should also note that for we are 
forced to account for the running of the anomalous dimension, rather than simply fitting 
to the free field dimension as suffices for mesons created by conserved currents. 

We also performed an analysis of the interaction terms allowed by the discrete and 
gauge symmetries of the model, and classified the lowest dimension and leading order in 
1/Nc terms that mix different states with the quantum numbers of the p,uj mesons in 
the presence of chiral symmetry breaking. In principle this model can now make new 
predictions for a number of physical quantities including 

• the mass spectrum and decay constants for the /ii,6i mesons, 

• the rates for hi — t- p + vr and bi u + tt as well as the D/S amplitude ratio for these 
decays, 

• the mass spectrum and vector and tensor decay constants for the low-lying excited p 
and oj meson states. 

Unfortunately this analysis requires rather complicated numerical analysis and as a result 
we postpone presentation of these results to a later paper [47]. 

Finally, we reiterate that although we have focused on the hard wall model, our con- 
siderations apply more generally to a broad class of dual models of QCD. The soft wall 
model and other models that modify the dilaton and/or metric require antisymmetric ten- 
sor fields of the type described here in order to incorporate hi /hi mesons. It would also 
be interesting to see whether the action described here for these fields can be derived for 
top-down models such as the Sakai-Sugimoto model using the techniques of open string 
field theory. 
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A Projection operators: definitions and useful relations 

We will often find it useful to consider the (anti) self-dual or longitudinal/transverse pieces 
of the two- form b^y. We can define these in terms of the operators and pH'-*-: 
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2 \ 2 
Useful relations among the projectors include 

p±^±p± ^ ^p± ^ p±^^pT = _ pll)pT _ (A.2) 

The projectors are idempotent, = P, while the projectors "pll'-*- have a non- 
standard normalization: = k'^V. The reason for this unaesthetic convention is that 
we do not wish projectors to harbor massless poles. To elucidate some of the potential 
subtleties, consider the derivation of the expression for n^'"^-*- in (2.13), starting with the 
definitions (2.2), (2.10), (2.11). By assumption, the states \b^^'^'^{k)) form a complete 
(properly normalized) basis for one-particle states that can be created by O^u from the 
vacuum. Therefore via insertion of the identity operator. 



(A.3) 



n,e 



where En{p) = yp^ + (m^"'')^. The sum over (on-shell) polarizations yields 



1 \^ -(n)a R (n)5 -y / (n)s2 i \ 

2 2^ eMt^a/3e(b) V'^ ^P^5f£\b) r = ) VplpVaju - {VplpPa]Pu - VulpPa]Pl,) 

e 

= {V^^^^^)p.,p. . (A.4) 
Using Cauchy's theorem, we can convert (A.3) to an integral over four-momentum, 



r ^4 1 /•("^)2 



p2 - (mj^)0^ 



Plugging this into (2.10) and using (2.11), we straightforwardly recover the expression for 
llT,T± (2.13). 

In the second step of (A.5) we used the fact that <f dzf{a)/{z — a) = ^ dzf{z)/{z — a), 
which holds provided / is a holomorphic function in the region bounded by the contour. 
This allows us to take the projector for each mode off shell, so that it is a common factor 
which may be pulled out in front of the summand. This step would have failed if V was 
a properly normalized projector, since then it would have a pole inside the contour of 
integration. 
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B General solution to the free equation of motion for Bmn 



The free equations of motion (3.8) for the two-form are 

ze^/" {-2ikpb^, + d,bp^) - 2ifibf,^ = , (B.l) 

ze^'P'^Kbp^ + 2fibp, = , (B.2) 

where we have used the Fourier-transformed fields 

bMN{k,z) = j d''xe'^-''bMN{x,z) . (B.3) 

Equation (B.2) yields the constraint 

h^.. = -^^e^'P'^k.bp, , (B.4) 

which, when plugged into the first equation of motion gives 



z 



' {k^kPbp, - k.kPbpp - k\,) + + fi\,. = , (B.5) 



or in terms of projectors, 

-z^V^ + ^z(P+ - p-)d, + /x^] bp^ = . (B.6) 

We can project (B.6) onto its self-dual and anti-self-dual parts. For the projections of 
the V'^b term, we write b = {P^ + P~)b and use (A. 2) as necessary. The resulting two 
equations may be rearranged to the form 

b^ = ^ (±^^5. + /x^ _ hh'^ {V^ - pll)6± . (B.7) 

We can thus plug the equation for 6~ into the equation for 6"*", deriving an equation for 6"*" 
alone, or vice versa. After some simplification and using {P-^ — pH)^ = /c^l, we arrive at 
the equations derived in [18] and quoted in (3.10): 

[z'^dl - zd, + fc^z^ - /i(^ ± 2)] = . (B.8) 

The general solutions are 

bpu{k,z) = Spy{k)zJ-p-i{kz) + Spy{k)zJ^,+i{kz) , 

bpp{k, z) = Ap^{k)zJ^p_+i{kz) + dp_^{k)zJ^_i{kz) , (B.9) 

where {S, s) and {A, d) are self-dual and anti-self-dual polarizations respectively, and J 
is a standard Bessel function. These solutions are appropriate (and convenient) for non- 
integer /i. For integer /i, the J-f^^i terms should be replaced by Yp_±i. We must remember, 
however, that we started with first order equations (B.7). These equations will be satisfied 
on the solutions (B.9) if and only if 

{A,d) = ^iP^-p\\){Srs) ^ iS,s) = ^{P^-p\\){A,d) . (B.IO) 
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Let us briefly comment on the bulk to boundary propagator for b. The form it takes 
depends on the sign of If > 0, then the S term of (B.9) dominates as z = e — )■ 0. 
The UV boundary condition matches 6+ onto a self-dual source, while the IR boundary 
condition sets b~ to zero: 



£M-l/2 



(B.ll) 



These conditions lead to the particular solution (3.25). On the other hand, if ;U < 0, the 
a term dominates. In this case h~ should match onto an anti-self-dual source at the UV 
boundary and h'^ should be zero at the IR boundary: 



= gi^i "m'^(^) ' ^M!^(^>^o) = , < 0) . 

These boundary conditions lead to the particular solution 



£\^'\B^{k,e) ' 
2 

{k^k^ CLpu ~ kj^k^app) 



Btik,z) 
■MB^{k,e) 



where 



(B.12) 



(B.13) 



B<{k,z) = zJ„|^|±i(/cz) - c^{k,zo)zJ\ 



with 



(kzo) 



(B.14) 



These expressions are appropriate for the non-integer /i case. For integer /i, the J_|^|±i 
should be replaced by i^^|=pi- 



C Dual sources, correlation functions, and matrix elements 

As usual in AdS/CFT, we determine the generating functional for the correlators of field 
theory operators by evaluating the supergravity action on the bulk-to-boundary propaga- 
tors of the dual fields. The UV boundary conditions on these propagators are determined 
in such a way that for a 4d source (j)o{x), the dual five-dimensional field has boundary 
condition^*^ (j){x,z = e) = (po{x). This determines the generating functional to be 

ZiMx)] = e*^™«™^'^^'^(-)='^o . (C.l) 

(Strictly speaking, this would be the leading saddle point approximation in a full quan- 
tum gravity dual description). We will usually work with momentum-space correlators. 
Functional differentiation in terms of momentum-space sources obeys 

^J—^J{p) = {27r)'5^'Hp-k) . (C.2) 

^"Up to powers of e, as discussed after equation (3.25). 
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We are mostly interested in correlators of the tensor operator qa^yq, but we often work 
with the pseudo-tensor or (anti) self-dual versions discussed in Section 2. Here we review 
the definitions of the sources for these operators and the relations between them. 

In the case /.f > the bulk to boundary propagator is determined by a self-dual source, 
P^S = S, which naturally couples to a self-dual operator in the field theory Lagrangian 
(in order to get a Lorentz scalar). The conjugate 5^,^ is anti-self-dual and naturally couples 
to the anti-self-dual, conjugate operator. These operators can be taken as 

0^,ix)=q{x)a,^-^^^q{x) , (C.3) 

and one indeed has (O^'^)* = 0^'~ . From them one can construct tensor and pseudo- 
tensor operators, 

Ol, = 0+ + = q{x)a^,q{x) , O^J = 0+ - O;, = q{x)a^,-i^q{x) . (C.4) 
These should couple to a tensor source and pseudo-tensor source, T^i, and V^j.u, given by 

V^,{k) = i {S^,{k) - S^,{-k)) S^,{k) = %,{-k) - V^.{-k) ' )■ ^ ■ ) 

Note that V,T satisfy V^u = ^e^u'^'^Tpa, as do the corresponding operators, O^^ and O^. 
The relative normalization between 7~, V and 5, S is fixed by the requirement that 

s^'^o^ + s^'^o;, = r'^'^oj, + v'^'^oj:^ . (C.6) 

When < on the other hand, the bulk to boundary propagator is given in terms of 
an anti-self-dual source a^,y. This source naturally couples to 0^'~, while a^,y couples to 
0^'~^ . The relationship of these sources to the tensor and pseudo-tensor sources introduced 
above is 



T^,uik) = ^(ttf^uik) + a^u{-k)) ^ apu{k) = Tpy{-k) -V^y{-k) 
V^uik) = \ (a^uik) - apu{-k)) a^u{k) = %u{k) + V^u{k) 



(^<0). (C.7) 



(Note that (C.5) and (C.7) do not imply a = S; these equations apply for different values 
of II, corresponding to different dual Lagrangians). 

We now derive the tensor-tensor two-point function from the free supergravity action. 
The same result was already derived in [18] , but we repeat it here for completeness. 

We begin by determining the generating functional for the two-point functions. First 
assume that /_f > 0. Evaluating the action on the bulk-to-boundary propagator (3.25), the 
bulk and IR boundary contributions to the action vanish, and we are left with the UV 
boundary term, 
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Plugging in (3.25), evaluated at z = e gives 



sd 



d^k B 



Agl J (27r)4e2^S+(A;,e) 
2^2/^-2 f d'^k B^{k,e) 



9l 



{2^^ye^^'B^{k,e) 



tr 



t Ci^A* UP Q 



+ c.c. 



(C.9) 



where we have used the fact that a self-dual tensor contracted with an ant-self-dual one 
vanishes. Since we are interested in the tensor-tensor two-point function, we choose to 
express the generating functional in terms of the tensor source. Using (C.5) with the 
constraint Vf^u = ^e^iu'^'^Tpa, a short calculation shows that 



h q^p-upq 



pp 



-T'^^i-k) [k^Tafsik) - 2 {ko^k^%p{k) - k^k"'Tya{k))] . (C.IO) 



Plugging this into (C.9), we arrive at (3.27): 



^sd 



d'^k B^{k,e) 



tiiT^^i-k) 



gl J {2TT)^k^e^PB+{k,e 
Now suppose that fi < 0. In this case the on-shell action is 

1 



r^^k)^ . (C.ll) 



5. 



sd 



d'^xti 



(C.12) 



Plugging in the bulk to boundary propagator (B.13), we find 



sd 



r B+{k,e) 



Agl J {2Tr)^e^M B~{k,e 
2£2|Mh2 |- ^4^ B+{k,e) 
{2TTYe^\^\B^{k,£) 



-tr< d' 



-al3 



+ c.c. 



9l 



tr 



(C.13) 



Using (C.7) with 

kpCi ^k^dpp 
and thus, 



- 2-f 

fiu — 2 P-'^ 'pot 



-T'^^ik) [k'Tapi-k) - 2 {k^k^T,p{-k) - kpk"'Tya{-k))] , (C.14) 



sd 



d^k 



B+{k,e) 



-tr 



[r^'^i-k) 



r^^(A:)} . (C.15) 



ff2 J (27r)4fc2£2|Hs-(A;,e) 
The two cases fj, > 0, fi < can be combined as follows. Define the functions 

B {k z) = \ ~ (^b{k,zo)zJ\^\_i{kz) , fj. non-integer, 

]^zY\p\^i{kz) - Cb{k,zo)zJ\p\_i{kz) , ^integer, 

B^(^k z) = I ^'^-Ia'I-i^^^) ~ Cb{k,zo)z.J\^\_^i{kz) , fi non-integer, 
)zYi^l+i{kz) - Cb{k,zo)zJipi+i{kz) , /i integer, 



(C.16) 
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with 



Cb{k,zo) 



(kzo) 



V|^|-l(A:2:o) 



fj, non-integer, 



/i integer, 



(C.17) 



for both positive and negative fi. Observe that when jj, > 0, Bn = BZ., and when /i < 0, 
Bn = -B<. Similarly, when > 0, B^ = Bt, and when fi < 0, B^ = B^. Thus, (C.ll) and 
(C.15) may be summarized as 



5sd = sgn(/i)- 



9l 



-tr <^ r 



■al3. 



-k) 



' al3,S^ ' a/3,57 



T'^k)}, 



{2tt)'^ k'^e^MBd{k,e) 

(C.18) 

valid for all fi. This is the expression we can functionally differentiate to compute the O^- 
two-point correlator. Due to the relation between V, T discussed above, the (O^^O^^) 

and (O'^O^^) correlators follow from the tensor-tensor correlator. 
The two-point correlator is computed via the usual prescription: 



{oUp)oU-'i)) 



sd 



Evaluating the right-hand side, restoring isospin indices, leads to 

Ae\^'\-^Bn{q,e 



(C.19) 



glq^e^MBd{q,e) 



fJ,U,pcr 



V, 



which gives a matrix element 



n 



T,T,ab/ 



^^.,pAk) = 5'^''sgn(/i) 
Given (2.11), we have H^'^H -- 



l£^M~'^Bn{k,e) 
ylk^e^MBdik,e) 

n^'^^, with 



V -V 



pp,pa 



pu^pa j ' 

(C.20) 



(C.21) 



n 



{k) = sgn(/i) 



glk^e^\i^\Bd{k,e) 



(C.22) 



Now, finally, we take the e — )• limit. We will first consider the case of non-integer /x, 
so that 

Bn{k,e) _ J-\f,\+i{ke) - Cb{k,zo)J\^\_i{ke) 



e^MBd{k,e) e2|H [j_|^^|_^(te)_Cfe(/c,zo)J|H+i(^^)] ' 



(C.23) 



For convenience, define Bn^d{k,z) = zB^^dik, z). We have that Bn/B^ = Bn/B^- Us- 
ing series expansions of the Bessel functions, we find that the leading e behavior in the 
denominator is 
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Meanwhile in the numerator we have 



2H-ir(|^ 

Hence the ratio takes the form 



^bik,zo) fclMl-i^lMhi(i + o(e2)) . (C.25) 



^(l^^l) 

Observe that all of the terms in the first series go as e to a negative or positive power; 
there are no terms that go as e° for non-integer The terms that diverge with e represent 
contact terms which we are not interested in. Hence, the first series can be ignored in the 
e — )■ limit. Meanwhile, only the first term of the second series survives the limit. Hence, 

lim /"^^'^^ ^ = contact terms - -^izMLk^Mcf,{k, zq) • (C.27) 



Next consider the integer fi case. We have 



(C.28) 



e^MBdik, e) e^l^l [Y\^\^^{ke) - c^ik, zo)J\^\+^{ke)] 
For integer n, the Bessel function K„ is defined by the series 

/ N 1 (n- j - 1)! /x\ 2i 2 , /x\ ^ , , 

^"(-)=-^(2) , M 2) +^1°H2)'^"^^)+ 

where ^/'o is the digamma function. For integer values, ijjo{n) = —'j + Hn-i = — 7 + X]j=i 
where 7 is the Euler-Mascheroni constant, and Hn a harmonic number. We note that the 
first series begins at order and terminates at order x^~^. The Jn term and the second 
series begin at order x". In the denominator then, the leading divergence is 

e^^i'^Bdik, e) = _ ^'^'^^l^^l' fc^lMl-igM^i (1 + 0(^2^^ ^ (^_3q) 

In the numerator, all terms in the first series of range from to el^l~^, and hence 

represent contact terms. The terms we are interested in, given (C.30), go as e^^^~^. They 
are the leading terms in J|^|_i, (both the explicit J\^\-i in (C.28) and the one contained 
in (C.29)), as well as the first term in the series on the second line of (C.29): 



Bn{k,e) = ■■■ + 



2 fke 
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where the • • • represent the contact terms. Thus we have that 



, Bn(k,e) 

hm —7—, ; = contact terms + 

IT 



22lH|/i|!(|/i| - 1)! 



2 / fee \ 1 

-log — - Cb{k, zq) {H\^\-i - 27) 

vr \ z / TT 



(C.32) 

In a standard renormahzation scheme with renormahzation scale Mr = this becomes 



BJk,e) 1 
..... ; = contact terms -\ — ^ 

6^0 e^MBd{k,e) 2^\> 

(C.33) 

In summary, (dropping the contact terms), we have 



hm ^1 1" — \ = contact terms -\ — — — - — ; — [vrcfef/c, zn) — log ik^f,^)] /c^'^' 



i^M}.c,(^, zo){klf\^\~^ , ^, non-integer. 



22|M|-2g2r(|^|) 
2^ra^g2|M|!(lMhl) 



[^"^(^' - log (^'^')] {klf^"^-^ , integer. 



(C.34) 

This reproduces (3.28) for /i > 0. Phigging in (C.17), we recover (3.31). 
D Spectrum of normalizable modes 

In this case we are interested in the normahzable part of the solution, (B.9). In the following 
we will consider the // > case only; an analogous treatment can be performed for the 
fi < case. For > 0, we restrict to the normalizable component of (B.9) by setting 
S = A = 0, resulting in 



4 , ■ 

s^u + '^k[fJ.^u]pk'^ 



zJ^^i^mz) . (D.l) 



Since bf,u = b'jlu + and b^z = -j^ef^^^'^kybpa, this implies 



bpy{k,z) = —s^„J^{7nz) + -^k^^Sy^pk^zJ^-i^mz) , (D.2) 



m 

1 2i 
baz{k,z) = e '''"'ki.Sp^zJJmz) = —Sauk''zJJmz) . (D.3) 

These solutions solve (B.6)-(B.8) with eigenvalue k'^ — t- m?. The IR Neumann-like bound- 
ary condition, b~{zo) = 0, fixes the eigenvalues to be m„ = x^_i^„/zo- For each eigenvalue, 
there is a corresponding self-dual polarization s^^{k), which contains six real degrees of 
freedom. 

We would like to evaluate (the quadratic part of) (3.7) on a sum over the eigenmodes. 
The boundary terms vanish since we have restricted to normalizable modes and imposed the 
IR boundary condition. To evaluate the bulk part of the action, we require the components 
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of db and -kb, evaluated on the solution (D.2). Working in momentum space we have 



{db)f_i,up = ^d[^fi^p] -i{kybup + Kbp^. + kpb^u) 



m 



(D.4) 



{db)^uz i\kfj^bi/z - kybpz) + dz 

4 2// 4 

= —k[pS„]pk''zJ^,{mz) + —SpydzJp{mz) + -^kipS^^pk^dzizJ^-^iimz)] 



2/i 

—SpudzJpimz) + Ak[p_Si,^pkP 



z 1 

— Jp{mz) + —Tdz[zJp-i{mz)] 



(D.5) 



while the components of ^ * ^ are 



2z 

ifj, 



2ifi 



ka^pupcr^ ^"^ S R^Ja(jnz) — {kpSyp + ki,Spn + kpSij,u)Jp{'m-z) 



pvz 



SpuJp{mz) 



pa 



SpuJaijnz) + -r-^k^kaEnyptje s fj^J n^\{rnz) 
mz fc^ 

2//2 

Ja{mz) + 2fiJp_i{mz) 

mz 



k°'Jp,-i{mz) 



4)U 

^f^i^ ~^ ~j^k\^pSy^pkP J p—i{rnz) . 



Therefore we see that 



while 



db-'-^kb 



, 



pup 



(D.6) 



(D.7) 
(D.8) 



db-'4*b 



2fi 



pvz 



—dzJJmz) + —Jp{mz) - J„_i(mz) 
m mz 



Spv~\~ 



+ 



k^Z 



m 



Jp{mz) + dz[zJp^i{mz)] - /iJ^„i(mz) 



^2 k[pSi,]pkP 



= k[pS^]pkPzJ^{mz) 

after making use of Bessel function identities. 
Then we have 

1 

2!3! 
1^ 



(D.9) 
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In the third step we plugged in some other, primed eigenmode for b and noted that the 
longitudinal part of b does not contribute to the contraction. We must multiply this result 
by —i when plugging into the action, and therefore the conjugate term gives a symmetric 
contribution, interchanging primed and unprimed modes. The IR boundary condition leads 
to an orthogonal spectrum, 

^ 

and therefore, defining canonically normalized polarizations 



gbrun 

the free action in the 6-sector takes the form 
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(fc-.-W) i^-^ {k^s(->p) . (D.13) 



Note that, since [b] = 3/2, we have [s] =5/2 and therefore [s] = 1. 

It is natural that only the longitudinal components of s should enter into (D.13): s is 
a self-dual tensor, so the transverse components are not independent degrees of freedom. 
The six real degrees of freedom contained in s may be taken, in the rest frame, as Re(soi) 
and Im(soi). According to our charge and parity assignments for b, discussed in section 4, 
the imaginary- longitudinal (equivalently, real-transverse) components of s represent hi/bi 
modes, while the real-longiudinal (equivalently, imaginary-transverse) components of s 
represent tensor oj/p modes, (2.3). Such modes can be given a canonical quadratic action 
by embedding them into s according to 



(D.14) 



As we discussed, at the level of the free action, these modes have identical spectra. 
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